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1. INTRODUCTION 

Using spectral imaging data acquired with the Airborne Visible Infrared Imaging 
Spectrometer (AVIRIS) from an ER-2 aircraft at 20 km altitude during various field 
programs, it was found that narrow channels near the center of the strong 1.38-pm water 
vapor band are very effective in detecting thin cirrus clouds (Gao et al., 1993). Based on 
this observation from AVIRIS data, Gao and Kaufman (1994) proposed to put a channel 
centered at 1 .375 (lm with a width of 30 nm on the Moderate Resolution Imaging 
Spectrometer (MODIS) (Salomonson et al., 1989) for remote sensing of cirrus clouds 
from space. The sensitivity of the 1 ,375-pm MODIS channel to detect thin cirrus clouds 
during the day time is expected to be one to two orders of magnitude better than the 
current infrared emission techniques. As a result, much larger fraction of the satellite data 
is expected to be identified as being covered by cirrus clouds; some of them so thin that 
their obscuration of the surface is very small. In order to make better studies of surface 
reflectance properties, thin cirrus effects must be removed from satellite images. 
Therefore, there is a need to study radiative properties of thin cirrus clouds, so that a 
strategy for correction or removal of the thin cirrus effects, similar to the correction of 
atmospheric aerosol effect (e.g., Kaufman and Sendra, 1988), can be formed. In this 
extended abstract, we describe an empirical approach for removing/correcting thin cirrus 
effects in AVIRIS images using channels near 1 .375 pm - one step beyond the detection 
of cirrus clouds using these channels. 

2. BACKGROUND 

Cirrus clouds consist of ice particles having different sizes and shapes. The 
"effective" particle sizes are usually greater than 5 pm. We illustrate the scattering and 
absorption properties of cirrus clouds through examples. Figure 1 shows two AVIRIS 
spectra (in reflectance units) acquired over areas covered by thick and thin cirrus clouds 
above Monterey Bay in California on September 4, 1992. The spectra below 0.5 pm are 
not shown here due to calibration problems. For each spectrum, the reflectances of ice 
particles in the 0.5 - 1.1 pm spectral region are constant with wavelength, because ice 
particle sizes are much greater than wavelengths. The larger apparent reflectances in the 
0.5 - 0.6 pm region are due to Rayleigh scattering. Beyond 1.1 pm, ice absorption 
effects, for example the ice bands centered near 1 .5 and 2.0 pm, are seen in the two 
spectra. The ice absorption effects near 1 .24 pm are weak. Near 1 .375 pm, both water 
vapor and ice particles have absorptions. The ice absorption effects near 1.375 pm are 
comparable to those near 1 .24 pm, because imaginary parts of ice refractive indices at 
both wavelengths are about the same. 

At present, accurate modeling of cirrus scattering and absorption effects is still 
very difficult. Our knowledge on cirrus particle size distributions and optical properties are 
not sufficient (King, 1993). Cirrus spatial inhomogeneity and adjacency effects are also 
difficult to be treated properly in radiative transfer models. In view of these difficulties, we 
propose to establish empirical relationships among remote sensing measures to 
characterize the main cirrus properties, and then to develop an empirical cirrus 
correction/removal technique. 
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Fig. 1 . Examples of AVIRIS spectra over thick and thin cirrus clouds. 


3. THE EMPIRICAL APPROACH 

For thin cirrus clouds, we assume that a homogeneous thin cirrus layer is located 
above a virtual "surface", which includes the effects of Rayleigh scattering and bottom 
surface reflection and scattering. With these assumptions, we have the following 
relationship (similar to descriptions of aerosol scattering, e.g., Fraser and Kaufman, 

1985): 

Px* = <>cX* (l) 

where p^* is the measured spectral radiance (normalized to reflectance units), p^ is 
the reflectance of the virtual "surface", p ^ is the atmospheric path radiance due to the 
cirrus cloud (radiation backscattered by the cloud), is the two way transmission (direct 
+ diffuse) through the cloud, and S ^ is the cloud scattering of upward radiation back to 
the surface. S - p_ « 1 for thin cirrus, and Eq. (1) can be simplified to: 

CA A 


= P cX + T cX P X 


( 2 ) 


In order to correct p^* for the cirrus effect, namely to derive p^ based on Eq. (2), cirrus 
path radiance p ^ and transmittance must be known. 

Analyzing AVIRIS data, we have found that cirrus path radiances p^ for 
AVIRIS channels between 0.4 and 1.1 pm are linearly related to P c ( 375 


p cA - P c 1.375pm ' ^ 


(3) 


where t^Q is less than 1 due to ice and water vapor absorption at 1.375 pm and can be 
derived from AVIRIS data themselves. For thin cirrus clouds, is usually greater than 

0.9. Correction of cirrus path radiances can be made relatively easily. Substitute Eq. (3) 
into Eq. (2), we obtain: 


T » Pi - P-i ■ P„ 1 aiSnm ^ *CO (4) 


l cX y X ~ ' H c 1.375pm 

Images of T ^ p^ are cirrus path radiance corrected images. Because of the large 

transmittance of thin cirrus, images of T ^ p^ are very similar to images of p^ 
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Fig. 2. A 0.648 pm AVIRIS image (top plot), 1.373 gm image (middle plot), and cirrus 
path radiance corrected 0.648 pm image (bottom plot). 

Corrections of the cirrus transmittances T ^ are practically more difficult than the 

corrections of path radiances. The solar radiation on the two-way path (Sun-cirrus-surface- 

cirrus-sensor ) encounters cirrus clouds at different spatial locations. Due to the non-local 

nature of T ^ , it is difficult to make accurate derivation of T « from AVIRIS images. 

CA C A 

Nevertheless, we are making progress in deriving a quantity from AVIRIS images that 

would characterize the main properties of T 

4. RESULTS AND DISCUSSION 

Figure 2 shows an example of our correction of thin cirrus effects. The top plot 
shows the 0.648-pm image. The image was acquired over Monterey, CA on September 
4, 1992. The image covers both land and ocean areas. A number of land surface features 
are blurred due to the presence of thin cirrus clouds. The middle plot of Fig. 2 shows the 
1.373-pm image over the same area. Only cirrus clouds are seen. Fig. 3 shows a scatter 

plot of p , „„„ versus p , ... for many pixels over the ocean. The straight line 
v r 1.373pm 1.043pm 
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Fig. 3. A scatter plot of apparent reflectances at 1.373 |im vs those at 1.043 jim. 



is a linear fit to all the points. The slope, 0.635, of the fitted line can be considered as 
good estimate of t qq (see Eq. 3). The bottom plot of Fig. 2 shows our path radiance 
corrected image. Most of the cirrus clouds over the ocean are removed. Land features are 
far less blurred than those in the top plot. 

As expected, the AVIRIS images in Figure 2 were contaminated by stratospheric 
aerosols resulted from Mt. Pinatubo volcano eruption in 1991. As a result, the minimum 
reflectance of the 1.373-pm is about 0.02, not zero. 

5. SUMMARY 

We have developed an empirical technique that is successful in removing the 
cirrus path radiances in AVIRIS images between 0.4 and 1.1 (im. The derivation of two- 
way cirrus transmittance factor is more difficult because of the non-local nature of the 
two-way transmittance. The authors would like to express appreciation to J. Von Den 
Bosch of Jet Propulsion Laboratory' for providing the AVIRIS data used in this study. 
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